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Asteroseismic modelling of the solar-like star (3 Hydri 
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Abstract We present the results of modelling the sub- 
giant star (3 Hydri using the seismic observational con- 
straints. We have computed several grids of stel- 
lar evolutionary tracks using Aarhus STellar Evolution 
Code (ASTEC, Christensen-Dalsgaard, 2008a), with 
and without helium diffusion and settling. For those 
models on each track that are located at the observa- 
tionally determined position of /3 Hydri in the HR di- 
agram, we have calculated the oscillation frequencies 
using Aarhus adiabatic pulsation package (ADIPLS, 
Christensen-Dalsgaard, 2008b). Applying the near- 
surface corrections to the calculated frequencies using 
the empirical law presented by Kjeldsen et al. (2008), 
we have compared the corrected model frequencies with 
the observed frequencies of the star. We show that af- 
ter correcting the frequencies for the near-surface ef- 
fects, we have a fairly good fit for both ^=0 and 1=2 
frequencies. We also have good agreement between the 
observed and calculated 1=1 mode frequencies although 
there is room for improvement in order to fit all the ob- 
served mixed modes simultaneously. 
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1 Introduction 

j3 Hydri is a G2-type subgiant star exhibiting solar like 
oscillations. It is often referred to as the future of the 
Sun due to having parameters (Table 1) close to those 
of the Sun, in addition to being at a later evolution- 
ary stage than the Sun. Being rather evolved, (3 Hydri 
exhibits mixed modes in its observed spectrum. This 
makes the star particularly interesting for asteroseismic 
studies as mixed modes carry more information about 
the core than do the regular p-mode frequencies. 



Table 1 Parameters of /3 Hydri as given in literature. The 
radius was derived from the interferometric angular diame- 
ter (North et al. 2007) and the revised Hipparcos parallax 
(van Leeuwen 2007). The mass was then derived from the 
mean density from asteroseismology. 



Parameter 


Value 


Reference 


AI/Mq 


1.085±0.028 


Kjeldsen et al. 2008 


r/r q 


1.809±0.015 


Kjeldsen et al. 2008 


L/L Q 


3.494±0.087 


Current work 


Teff(K) 


5872±44 


North et al. 2007 


[Fe/H] 


-0.08±0.04 


Santos et al. 2005 


Ai^o(/iHz) 


57.24±0.16 


Bedding et al. 2007 



An extensive analysis of this star was done by Fer- 
nandes & Monteiro (2003), who emphasized the possi- 
bility of employing seismic constraints to remove par- 
tially the parameter degeneracy that exists when only 
the non-seismic observational constraints (such as T e g 
and luminosity) are used. Those constraints alone are 
not enough, as one can obtain the same location of a 
stellar model in the Hertzsprung- Russell (HR) diagram 
using different stellar parameters (such as mass, X, Y, 
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etc.). Fernandas & Monteiro (2003) used the large 
frequency separation, Ai'o, obtained by Bedding et al. 
(2001), to show how to constrain the mass indepen- 
dently and also note the need for individual frequencies 
to further constrain the age of the star using the char- 
acteristics of the frequency spectrum that are related 
to the stellar core, such as mixed modes. 

The first preliminary comparison of the individual 
observed frequencies and model frequencies was done by 
Di Mauro et al. (2003). They presented a model with 
large and small frequency separation that are within 
the limits derived from the observations (Bedding et al. 
2001). They showed that the match between the model 
and observed frequencies are satisfactory except for the 
1=1 modes, some of which are affected by avoided cross- 
ings. This emphasizes again the importance of accurate 
asteroseismic observations and detailed analysis in or- 
der to evaluate the stellar interiors. 

Here, we use the latest asteroseismic observational 
constraints (Bedding et al. 2007) consisting of individ- 
ual frequencies including some modes which are identi- 
fied to be possible mixed modes. We present the meth- 
ods to search for a best model and the resulting best 
fits within the parameter space of our survey. 



2 Methods 

We started calculating the grids of evolutionary tracks 
with a wide range of parameters and large steps of in- 
crement. Analyzing the first grid we have selected the 
best models, around which we have computed denser 
grids. The initial parameters of the grids are given in 
Table 2. Diffusion and gravitational settling of helium 
are added in the third grid, in which not all the tracks 
have been successfully completed. We have carried on 
our analysis with the tracks that did not have conver- 
gence problems. 

On each track in the grids, we have selected the mod- 
els having parameters within the observational uncer- 
tainty limits. We have calculated the oscillation fre- 
quencies of those models and compared with the ob- 
servations. For the purpose of comparison we used 
the frequencies resulting from dual-site radial velocity 
observations with HARPS and UCLES spectrographs 
(Bedding et al. 2007). 

Before comparing the calculated frequencies with the 
observed ones we have applied near-surface corrections 
to the calculated frequencies. The correction is needed 
due to the fact that existing stellar models fail to rep- 
resent properly the near-surface layers of the solar-like 
stars, where the turbulent convection takes place. This 
affects the high frequencies most; thus the correction 



should be applied in a way that low frequencies are 
much less affected. 

The situation is the same for the Sun, and the dif- 
ference between observed and calculated frequencies is 
shown to be well approximated by the empirical power 
law given by Kjeldsen et al. (2008) as 



v bs{n) - fbest(«) = a 



(1) 



where v ^ s are the observed 1=0 frequencies with radial 
order n, fb os t are the corresponding calculated frequen- 
cies of the best model, and Vq is a constant frequency 
chosen to be the frequency corresponding to the peak 
power in the spectrum, which is taken as 1000/iHz for 
(3 Hydri. Kjeldsen et al. (2008) used the solar data 
and models to calibrate the exponent b, which is cal- 
culated as 4.90 for the Sun. Using this solar b value, 
and calculating a for each model, we have applied the 
near-surface corrections to the model frequencies. The 
corrections for mixed modes are probably less than for 
pure p modes, as discussed by Kjeldsen et al. (2008), 
and so we have set them to zero. We have then selected 
the best models performing a chi-squarc minimization 
test for the goodness of the fit as 



.model 



(n) 



,obs 



" ) 



er(i/° bs (n)) 



where N is the total number of modes included, i/f hs 



(2) 



(n), 



and i/ ; model (n) are the observed frequencies, and the cor- 
rected model frequencies, respectively, for each spher- 
ical degree I and the radial order n, and a represents 
the uncertainty in the observed frequencies. 



3 Results 

Properties of our best models are given in Table 3, 
with the corresponding so-called echelle diagrams in 
Figures 1 and 2 (see Christensen-Dalsgaard (2004) and 
references therein for the explanation of an echelle di- 
agram). The increase in the systematic difference be- 
tween the observed and model frequencies with increas- 
ing frequency can be seen in the left panels of the echelle 
diagrams. 

The echelle diagrams include all the observed fre- 
quencies, while those having identified modes are con- 
nected by a solid line for clarity. The observed frequen- 
cies which are not connected by a line and shown with 
smaller symbols are the ones that are not assigned a 
mode by the observers (Bedding et al. 2007); however, 
they note that those frequencies will include some gen- 
uine modes although some of them might be sidelobcs 
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Table 2 Parameters used to compute the evolutionary tracks 



Parameter 


Grid 1 


Grid 2 


Grid 3 


M/Mq 


1.04-1.10 (with steps of 0.01) 


1.076-1.084 (with steps of 0.002) 


1.076-1.084 (with steps of 0.002) 


z/x 


0.018-0.022 (with steps of 0.001) 


0.018-0.022 (with steps of 0.001) 


0.018-0.022 (with steps of 0.001) 


Y 


0.23, 0.27, 0.28, 0.30 


0.276-0.284 (with steps of 0.002) 


0.276-0.284 (with steps of 0.002) 


Mixing length 








parameter (a) 


1.4-1.8 (with steps of 0.2) 


1.75-1.85 (with steps of 0.025) 


1.75-1.85 (with steps of 0.025) 


Diffusion & 








gravitational settling 


None 


None 


He 



Table 3 Parameters of the best models 



Parameter 


From grid 2 


From grid 3 


M/Mq 


1.082 


1.082 


z 


0.01346 


0.01266 


Y 


0.278 


0.284 


a 


1.825 


1.775 


Age (Gyr) 


6.447 


5.712 


R/Rq 


1.814 


1.806 


L/L Q 


3.444 


3.473 


T cf r(K) 


5844 


5869 


a 


-2.76 


-3.76 


x 2 


7.71 


11.93 



before correction 



after correction 
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Fig. 1 Echelle diagram of the best model without gravita- 
tional settling or diffusion. Left (right) panel shows the case 
before (after) applying near-surface corrections. Stars de- 
note the observations, squares correspond to the model fre- 
quencies with 1=0, triangles to the model frequencies with 
1=1, and diamonds to the model frequencies with 1=2. Dot- 
ted vertical lines indicate the value of the large frequency 
separation, Ai/q 



or noise peaks. The size of the largest symbols cover 
approximately 1.5-c uncertainty to both sides on the 
horizontal axis. Three of the identified 1=1 frequencies 
that fall to the left side of the 1=1 ridge are identified 
as mixed modes by the observers (Bedding et al. 2007). 
The model without He settling or diffusion reproduces 
the observed frequencies very well except for the lowest 



mixed mode. It is seen from Fig.l that some of the 
unidentified observed frequencies also match the model 
frequencies quite well. When the diffusion and settling 



before correction 



after correction 
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Fig. 2 Echelle diagram of the best model with He settling 
and diffusion. The symbols are used in the same way as in 
Fig. 1 

of helium is added, the agreement between the model 
and the observations is less strong than the previous 
case. Although the lowest mixed mode is reproduced in 
this model, the mixed mode having the sharpest charac- 
ter (largest departure from the 1=1 ridge) is not. The 
highest mixed mode has not been reproduced by any 
model frequency with 1=1, suggesting that 1=3 frequen- 
cies might be investigated to search for a match to that 
frequency The positions of both models on the sub- 
giant branch are shown in Fig. 3. The model with 
diffusion is younger due to having completed the main 
sequence phase faster, as the hydrogen mass fraction 
decreases faster owing to diffusion and settling bring- 
ing helium to the core. Furthermore, the hook shape 
in the evolutionary track with diffusion is due to the 
fact that the star had grown a convective core up to 
an extent where it contained almost 4% of the stellar 
mass. 
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Fig. 3 Evolutionary tracks with He settling and diffusion 
(dashed curve, the best model marked with a circle) and 
without (solid curve, the best model marked with a cross) 
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4 Conclusions 



Our results justify that the empirical power law repre- 
senting the effect of near-surface layers in the Sun works 
for (3 Hydri as well. Our best models with and without 
He settling and diffusion reproduce the observed 1=0 
and 1=2 modes well; however, the fit at 1=1 modes is 
relatively poor due to the mixed modes, but still satis- 
factory. It is important to investigate the possibility of 
any of the observed modes being an ^=3 mode. 

Furthermore, our results are in agreement with the 
findings of Fernandes & Monteiro (2003), who derived, 
through the HR diagram analysis, the mass to be 
l.lOlg'Q^ -^O) the helium abundance to be between 
0.25 and 0.30, and the stellar age to be between 6.4 
and 7.1 Gyr. We, however, note that the age of our 
model with helium diffusion is less than the lower limit 
of the cited result. Employing Av , they found the mass 
to be 1.09±0.22 M©, noting that the large uncertainty 
is to be improved with the improved accuracy of the 
observations. Further analysis may be carried out to 
investigate the effect of convective core overshooting as 
our models had convective cores at some earlier stage 
in their evolution. 
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